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ABSTRACT 



We investigate the effect of the environment on the Faber Jackson (FJ) rela- 
tion, using a sample of 384 nearby elliptical galaxies and estimating objectively 
their environment on the typical scale of galaxy clusters. We show that the in- 
trinsic scatter of the FJ is significantly reduced when ellipticals in high density 
environments are compared to ellipticals in low density ones. This result, which 
holds on a limited range of overdensities, is likely to provide an important obser- 
vational link between scaling relations and formation mechanisms in galaxies. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: fundamental 
parameters 
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Introduction 



The Faber Jackson (FJ) is the first scaling relation discovered for elliptical galaxies. 
Already described by Morgan & Mayall (1957), although in a qualitative way, it was 
given its first quantified form 20 years later by Faber & Jackson (1976) who, on the basis 
of a handful of nearby early type galaxies, were able to proof the existence of a power 
law relation linking luminosity (Lb) to central velocity dispersion (00). Soon thereafter 
Kormendy (1977) found a second scaling relation which holds for elliptical galaxies, relating 
the effective surface brightness (/x e ) to the effective radius (Re). The Kormendy relation 
was refined ten years later for ellipticals by Hamabe & Kromendy (1987) and in that same 
year Dressier et al (1987) and Djorgovsky & Davis (1987) discovered a more general relation 
(the Fundamental Plane, FP) linking Log R c , Log er and // e . 

The scaling relations are powerful tools that can be used to derive galaxy distances 
and, even more important, constitue an invaluable observational bench mark for theoretical 
models. It is especially for this latter reason that they have been the subject of much 
interest since their discovery. Understanding origin and nature of the scaling relations is a 
fundamental quest for any successful theory of galaxy formation which is expected to be 
able to predict the observed slope, scatter, possible variation (as a function of luminosity, 
wavelength, environment) and evolution (with z). The much narrower scatter displayed by 
the FP with respect to the FJ and Kormendy relations made the former one to become 
rapidly more attractive than the latter two which were easily interprete d to be partia l 



repres e ntations (projections ) of t 



( 119871 ): iDjorgovski fc David ( 119870 



r e FP onto a l ower d imensional space 



Faber et al. 



Dressier et al. 



(jl987[ ): Ide Zeeuw fc Franxl (jl99ll ). 



Very recently the FJ appears to have captured again attention on the theoretical 
point of view, as Sanders (2010) has claimed it to be more fundamental and universal 
than the FP within the context of MOND (modified Newtonian dynamics). This finding 
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is expected to motivate renewed interest in the FJ, which so far has not been largely 
investigated. There is no much work which has been carried out on the FJ relation 



if one excepts stu dies which ha v e provided evidence 



at low 



(2005); 



uminosity 



Tonrvl ( 



198 



Matkovic fc Guzman (2005); 



LaueretaL 



(120071 ) 



Davies et al 



(119831) : 



Ber lardi ct a 



Von der Linden et al. 



(12007 



or a decre a se of it s stepness 



Held et al. 



(2006); 



(1992); 



Fritz et al. 



Desroches et al 

— n — i — 



( 120071 ): 



Kourkchi et al.l (120121 ) and studie s 



devotee 



( 120051 ): 



to investigate 



Bernardi et al. 



( 120061 ): 



Desroches et al. 



(120071 ): 



and redshift on it 


Fri 


tz et al. 




Nigoche-Netro et al. 


( 


2010. 


2011) 



At variance with the FP for which the effect of the environment has been largely 



investigated, although with ra t her conflicting results 



Mar quez fc Moled ( 



(2003); 



20051 ): 



1996 



Evistingeeva al. 



Denicolo et al 



4 studies exist 



Zieglerl (12005 



1999); 



(12002 ): 



(2005); 



de la Rosa et al 



de Carvalho fc Diorgovskil (119 92); 



(2001); 



Treu et al. 



Gonzales-Garcia fc van Albada (2003); 



D' Onofrio et al 



); 



Fritz et al 



(12005 



(2008) 



20091 ): 



(2001); 



Bernardi et al. 



-teda et al. 



La Barbera et a 



Fritz &; Zieglerl (120091 ) which have 



( 2004 



(l2O10h . so far only 



looked for possible effects induced by the environment on the FJ relation of rather distant 
(z e [0.2 - 0.7]) early-type galaxies, without finding however any strong evidence for them. 

According to the standard cosmological paradigm, structures in the present day 
Universe have formed through a hiearchical scenario process predicting rather different 
assembling time scales and evolutionary paths for galaxies in high and low density regions 
(Baugh et al. 1996; Kauffmann & Chariot 1998; Somerville & Primack 1999; Kauffmann et 
al. 2004). Environment is thus expected to play a relevant role in shaping galaxy properties 
and is likely to leave its imprint in the scaling relations as well. This is the reason which has 
motivated the above mentioned studies (mostly concentrated on the FP) and the present 
work devoted to investigate the effect of the environment on the FJ relation, using a sample 
of 384 nearby ellipticals and estimating their environment on the typical scale of galaxy 
clusters. 
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The structure of the paper is the following: in §2 we present the sample, in section 
§3 we derive the FJ relation for the whole sample and for its bright and faint components 
and test the robustness of our results accounting both for errors on cr and m#, in §4 we 
illustrate the method that we have used to estimate the environment, in §5 we show that 
the scatter of the FJ gets largely reduced in high density environments and increased in low 
density ones and that this difference is neither induced by errors on <jq nor by luminosity 
difference between the samples, in §6 we show that the scatter of the FJ relation increases 
with decreasing density in overdense environments and decreases with increasing density in 
underdense environments, in §7 we draw the conclusions. 

A Hubble constant of H = 70 km s" 1 Mpc -1 is adopted throughout. 



2. The sample 



We have extracted the sample of elliptical galaxies from HyperLeda iPaturel et al 



(120031 ) requiring them to have a measured, m B < 15.5, v T G [2000 - 10000] km s -1 , and 
5 > ° . The constraints have been imposed to homogenize a sample that otherwise would 
not satisfy any predefined selection criteria, having been drawn from a compilation of 
available data. The limit in m-Q makes the sample to become flux limited, the lower limit in 
v r reduces distance uncertainty due to the contribution of peculiar motions, the upper limit 
in v r keeps contained the increase of luminosity with increasing distance and the limit in 
S ensures homogeneous sampling of the environment for all ellipticals in the sample, since 
the catalog that we have used to estimate the environmnent (see next session) is limited to 
5 > -2.5 °. 

The selection based on the above described prescriptions produced a sample of 
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384 elliptical galaxies, that are listed in Table 1 Q. For each elliptical, Table 1 reports 
identificator (column 1), equatorial coordinates RA J2 ooo an d -Decj 2 ooo (columns 2 and 
3), central velocity dispersion o" and related uncertainty Acr (columns 4 and 5), total 
apparent B magnitude tub and related uncertainty Atub (columns 6 and 7), radial velocity 
v r corrected for the Virgo infall flow (column 8), number of neighbours A^ neig h detected 
within the typical group/cluster scale (the details on the neighbour search method will be 
given in the next section). 

Figure 1 shows the distribution of our sample ellipticals (filled circles) in the v T , Mb 
plane together with the curve corresponding to the faintest observable Mb in a sample 
limited to me = 15.5 (which is the limit that we have imposed to our sample). The curve 
allows one to visualize the well known effect induced by distance on luminosity in a flux 
limited sample: since the minimum observable luminosity increases with distance, the 
farthest galaxies will be, on average, also the brightest ones. From Fig. 1 we see, however, 
that the distance effect constrains only the minimum observable luminosity and that its 
real entity depends quite strongly on the galaxy distribution in the v T , Mb plane, which will 
never be completely uniform, due the intrinsic dumpiness of galaxy distribution (i.e. the 
presence of clusters and groups) and in the present case also to the possible lack of some 
data which we cannot exclude as we are dealing with a sample drawn from a compilation of 
available data. 

The effect of distance on the luminosity can be perceived better from Fig. 2 which 
shows the absolute magnitude (Mb) distribution of ellipticals in the whole sample (dotted 
histogram) overimposed to the Mb distribution of ellipticals in 4 bins of increasing radial 
velocity, each spanning a 2000 km s _1 range. Definition and characteristics of the bins 
are given in Table 2 providing for each bin (column 1), radial velocity range (column 2), 



online only table. 
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total number of ellipticals (column 3), mean (column 4) and median (column 5) value of 
the Mb distribution (these latter to be compared with -20.80 and -20.87, which are the 
corresponding values of the whole sample). Figure 2 evidences rather clearly the progressive 
shift of the Mb distribution of ellipticals as a function of the increasing distance and from 
columns 4 and 5 of Table 2 we see that this distance effect induces an average increase of 
1.2 magnitudes between the farthest and the nearest ellipticals in our sample. 



3. The FJ relation for the whole sample 

The FJ relation for the total sample is shown in Fig.3. Superimposed to the data 
(open circles) is the best fit (green continuos line) derived taking into account errors on 
do (displayed in Fig. 3 as well) and corresponding to Lb oc <Jq°. The red dashed (Lb oc 
<Tq' 6 ) and blue dotted (Lb oc <Jq' 2 ) lines instead represent the best (weighted for errors on 
(T ) fits for the 206 bright (Log(L B /L B0 ) > 10.5) and the 178 faint (Log(L B /L B0 ) < 10.5) 
ellipticals in the sample, with Log(LB/LB©) = 10.5 being the value below which data in 
Fig. 3 start do deviate progressively from the fit and to extend towards low values of do- 

From Fig. 3 we see that the dispersion of the data around the fits is quite large with the 
average scatter (RMS), cr(o"o), amounting to 51.1 ± 1.2 km s _1 for the whole sample and to 
49.3 ± 1.8 km s _1 and 52.3 ± 1.7 km s _1 respectively for the bright and faint subsamples. 



The less steep increase of Lb with cr displayed b y faint 



confirms previous evidence ob t ained by several authors 



Held et al. 



(1992 



Desroches et al. 



Fritz et al. 



( 120070 : 



(12005 ) 



Lauer et al 



Tonrv 



Matkovic k Guzman (2005 



( 120070 : 



ellipti c als in our samp le 



198 



Von der Linden et al. 



); 



Davics et al. 



Bernardi et al. 



(2007 



1983); 



(2006): 



Kourkchi et al. 



( 120121 ) on different samples which find theoretical justification iDekel fc Silk! ( 119861 ) in the 



expected link beween the FJ relation slope and the amount of dark matter in elliptical 
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galaxies. 

Some caution should be given however to slope values derived by means of fits in which 
errors on <r (column 5 of Table 1) have been taken into account, as they might result 
somewhat artificially increased if data with large o"o errors (which weight less) are mostly 
found below the fit lines. This seems actually to be the case in our sample (cfr. Fig. 3) and 
from Fig. 4 we can visualize better the effect induced by the weighted fit on the slopes. 
The distributions of the relative o"o errors (Afr /o"o) of data laying below (continuos) and 
above (dotted) the best fit lines for the unweighted (left) and weighted (right) fit of the 
whole sample (upper panels) and of the bright and faint subsamples (middle and lower 
panels) presented in Fig. 4 show that the effect produced by the weighted fit is to reduce 
the numerical dominance of ellipticals with small A<t /(t above the fit lines increasing, as a 
consequence, the slope of the FJ relation. Had we derived the slopes with an unweighted fit 
we would have obtained L B oc a^ A for the whole sample, L B oc o"q 7 for the bright and L B 
oc dp' 7 for the faint subsample (with average scatters cr(a ) respectively equal to 51.9 ± 1.2 
km s -1 , 49.5 ± 1.8 km s" 1 and 52.9 ± 1.7 km s' 1 ). 

It is worthwhile to stress that a weighted fit will always produce an artificial steepening 
of the slope (whatever the size of the errors) if data with large errors are preferentially found 
below the fit. To prove that we have derived the FJ relation for the 314 ellipticals having 
Acro/oo < 0.1, finding that the weighted fit induces an increases of the slope from 3.8 to 4.1 
for the whole sample, from 5.4 to 5.8 and from 3.0 to 3.3, respectively for the bright and 
faint subsample. The average dispersion of the data diminishes a little and sets around 48 
km s _1 , the exact value depending on the kind of fit and sample. The slight reduction of 
cr(cxo) is expected as we have excluded data with large relative errors on <t (he. Aa /ao > 
0.1) which being less accurate are more likely to deviate more strongly from the fit. 

Table 3 summarizes all the results described above. In column 1 we list the sample kind 
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(whole, bright or faint), in column 2 the number of ellipticals in each sample, in column 3 
the kind of fit (either unweighted or weighted), in column 4 the FJ relation slope (a) with 
related uncertainty, in column 5 the average scatter of data around the best fit line (er(er )) 
with related uncertainty. Figures in Table 3 do not allow us to establish the exact value for 
the slope of the FJ either of the whole sample or of its bright and faint components, but 
allow us to confirm the presence of two distinct (luminosity dependent) components in the 
FJ relation characterized by a slope which is steeper for bright (Log(L B /L B0 ) > 10.5) than 
for faint (Log(L B /L B0 ) < 10-5) elliptical, as the difference between the slopes holds (and 
is larger than the errors) whatever the kind of sample (either whole or whole with small 
Actq/o-q) and of fit (either unweighted or weighted). 

Finally, to check the effect of possible errors on m B on the derived FJ relation, we have 
randomly added or subtracted to each m B (column 6 in Table 1) either its real (column 7 in 
Table 1) or average error (computed on the whole sample). We have repeated this operation 
300 times for both cases, thus obtaining 2 sets of data each including 300 simulated 
samples. We have subsequentely derived the FJ relation (weighted fit) for each simulated 
sample in each set and grouped the results together, to obtain the total distribution of 
er(ero), a, and a max - a m i n (i.e. twice the maximum uncertainty on the slope a). These 
distributions (normalized to the total number of simulated samples) are shown in Fig. 5, 
where continuos and dotted histograms refer respectively to samples obtained by random 
addition or subtraction of the average or of the real error on m B . Upper panels refer to 
the whole simulated samples, middle and lower panels to the bright and faint subsamples, 
extracted from the previous ones. Since the artificial random increase or decrease of m B 
implies a decrease or increase in L B , the number of ellipticals in the high or low luminosity 
subsamples is not anymore constant but varies as a consequence of the variation in L B . We 
find a median value of 204 and 180 ellipticals and of 210 and 174 ellipticals (to be compared 
with the real value of 206 and 178) for the high and low luminosity subsamples when the 
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average or real error on tub is respectively randomly added or subtracted. 

The arrow placed on each distribution of Fig. 5 indicates the corresponding value 
obtained in the real case (i.e. the weighted fit for the whole sample or for the bright or for 
the faint subsample, cfr. figures in Table 3, lines 4,5 and 6) and allow us to state that the 
effect of the error on m B is an increase in the average dispersion (a(a )), a steepening of 
the slope (a) and an increase of the uncertainity of this latter quantity (a max - ctmin)- From 
Fig. 5 we see that the effect of the error on me is stronger for the whole sample than for 
the faint and bright subsamples, which is not surprising as we have shown (in sect. 3) that 
the FJ relation of the total sample can be interpreted as due to the combination of two 
distinguished (luminosity dependent) components. Variations in L B of each elliptical in the 
sample are thus expected to produce a stronger (amplified) effect on the FJ relation of the 
whole sample than on the FJ relations of the separate (luminosity dependent) components. 
From Fig. 5 we also see that the effect of the error on m B is stronger for the faint (lower 
panels) than for the bright (middle panels) subsample. This is not surprising too, since 
variations in Lb are expected to influence more strongly fits which have a less steep slope. 

In Table 4 we list for each sample (column 1), the error Am B , (either the average value 
for the whole sample or the range within which the true value is found) which has been 
randomly added or subtracted to tub (column 2), the mean value of the average scatter of 
the data around the fit <a(a )> and its RMS (columns 3 and 4), the mean value of the 
slope <a> and its RMS (columns 5 and 6), the mean value of a max - a min (here indicated 
as <Aa>) and its RMS (columns 7 and 8). Comparing the mean values listed in Table 4 
(columns 3,5 and 7) with values obtained for the real sample and subsamples (arrows in 
Fig. 5 corresponding to figures listed in Table 3, lines 4,5 and 6) we see that on average the 
effect of the error on itlb can be considered moderate. Data in Table 3 and Fig. 5 allow 
us to conclude that even if the tub of each elliptical would vary of a quantity respectively 
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equal to the average or real error, which is surely an overestimate of what may happen in 
the real case, we could however confirm the value of the average dispersion of the data 
(a(ao)) around the best fit which would get only slightly increased. The presence of two 
distinct (luminosity dependent) components would be confirmed as well, since both slopes 
would get somewhat increased but remain still well distinguished (their difference being 
larger than their errors). 



4. The environmnent 



To evaluate the environment of each elliptical galaxy in our sample we have applied the 
neighbour search code of Focardi & Kelm (2002) to the Updated Zwicky Catalog (UZC) 



Falco et al. 



(1999) 



galaxies that cove rs the entire northern sky 



Falco et al. 



fll999f ) to be 96% complete for 



UZC is a wide angle 3D catalog of nearby 
down to a declination of -2.5°, and is claimed 
galaxies brighter than m-Q = 15.5. 



The neighbour search code is a versatile tool which can be applied to 3D catalogs either 
to produce galaxy samples characterized by different environment or to estimate galaxy 
environment on different scales and depth. 

Detailed description of the code can be found in Focardi & Kelm (2002), together 
with the results of the first application of the code to UZC which has produced a large 
homogeneous sample of compact groups (UZC-CGs, Focardi & Kelm 2002). The code has 
been subsequentely applied to UZC to produce a sample of bright isolated galaxy pairs 
(UZC-BPGs, Focard i et al. 2006) and a small sample of very isolated bright elli pticals 



Memola et al. 



(1200a ) . It has also been applied to the 2dFGRS 



Colles et al. 



(120011 ) in order 



to perform a detailed analsys on the luminosity /environment /spectral type relation for 
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galaxies iKelm et al.l ( 120051 ) 



When used simply to detect neighbours, as in the present case, the code needs only 
two input parameters which are the maximum projected distance (Ai2) and radial velocity 
difference (|Av r |) between each elliptical in the sample and its possible neighbours (from 
UZC) and includes obviously, a cross check on coordinates, v x and m-Q to avoid spurious 
detection of the elliptical as possible neighbour of itself. 

We have set AR = 1.5 Mpc and \Av T \ = 1000 km s _1 , a choice which has allowed us to 
estimate the environment on the typical scale of galaxy clusters. 

Figure 6 shows the distribution of the number of neighbours (iVneigh) when the whole 
sample is divided into the 4 bins of increasing radial velocity, that we have defined in §2 
and whose characteristics are summarized in Table 2. Distributions appear rather different, 
which is not unexpected as we have already shown (cfr. Figs. 1 and 2, Table 2 columns 
4 and 5) how distance effect produces an increase of luminosity with increasing distance. 
Thus, as the number of galaxies decreases with increasing luminosity we would expect, 
on average, a decreasing number of neighbours with increasing distance. Figure 6 shows 
that this is not exactly the the expected decrease in iV nei g h is evident only in the 

fourth bin, while both bins II and III show an anomalous tail (large values of iV neig h), 
which is due to the presence of several galaxy clusters and groups (belonging respectively 
the Perseus-Pisces and Coma supercluster) of which some ellipticals in bins II and III are 
members. 



5. The effect of the environment on the FJ relation 

To look for possible effects induced by the environment on the FJ relation one must 
compare ellipticals in high and low density environments, with density being as large and 
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as small as possible but leaving however a number of ellipticals in each environment that is 
not too small. 

An objective way to define these extreme environments can be obtained relating the 
number of neighbours requested to enter each sample to the median value of the A" ne igh 
distribution, this latter computed separately for each bin to account for effects related to 
distance and non uniformity in the galaxy distribution. 

We find that defining as high or low density an environment characterized by a number 
of neighbours equal or larger or equal or smaller than 5 times or 0.2 times the median value 
of iVneigh provides us with two subsamples of 26 and 36 ellipticals (respectively including 
7% and 9% of the whole sample). Since the median value of iV ne i g h is equal to 8 in first bin, 
to 9 in the second and third bin and to 4 in the fourth bin, our definition of high density 
environmnent implies iV ne igh > 40 (in bin I), A^ nei gh > 45 (in bins II and III) and iV ne igh > 20 
(in bin IV), while our definition of low density environment implies iV ne i g h < 1 (in bins I, II 
and III) and iV ncigh = (in bin IV). 

Figure 7 shows the FJ relation for ellipticals in the high (red filled circles) and low 
(blue open squares) density environments defined above. Each data is displayed with its 
own error (both in m B and in cr ), while the red continuos and blue dotted lines represent 
the best weighted fits respectively giving L B oc cr^ 79 and L B oc o-q 75 - It is evident that 
ellipticals in low density regions distribute in a much more dispersed way than ellipticals in 
high density ones. The o"(o" ) of the FJ in low density environments (67.0 km s _1 ) is in fact 
almost twice the one found in high density ones (33.7 km s -1 ) and even larger than what 
found for the whole sample (51.1 km s _1 ). 

Weighted and unweighted fit give exactly the same slope (a ~ 3.8) in high density 
environments, while in low density environments the slope is steeper in the weighted fit (a 
~ 3.8) than in the unweighted one (a ~ 3.3), due to the effect produced by the dominance 
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of data with larger o"o errors below the fit line (as discussed in sect. 3). 

One might then argue that the large er(<7 ) displayed by ellipticals in low density 
environments could be induced from their large o"o errors, but if we exclude from the sample 
the 6 ellipticals with the largest relative errors (Aa /a > 0.15) we still get large values 
for <t(<To) (67.9 km s -1 and 68.5 km s _1 respectively for the weighted and unweighted fit). 
Moreover, from Fig. 7 we see that the luminosity distribution of ellipticals in high and low 
density environments is rather similar, implying that the larger value of cr(a ) displayed by 
ellipticals in the latter sample cannot be attributed to a luminosity effect linking intrinsic 
scatter to luminosity, as already evidenced for the FP (see e.g. Bender et al. 1992; Hyde & 
Bernardi 2009). The fraction of faint (Log(L B /L B0 ) < 10.5) ellipticals is in fact almost the 
same in high (15/26, ~ 58 %) and low (20/36, ~ 56 %) density environments, implying that 
the larger o"(o"o) displayed by ellipticals in low density environmnent cannot be attributed 
to a larger content of low luminosity ellipticals. 

An accurate comparison of the luminosity distribution of ellipticals in high and low 
density environments can be obtained from Fig. 8, which shows the normalized M B 
distribution of ellipticals of the two samples. From Fig. 8 we see that ellipticals in 
high density environments (shaded histogram) cover a somewhat larger luminosity range, 
extending at both sides of the distribution while ellipticals in low density environments 
dominate at M B ~ —21. The median value of the distributions is exactly the same (M B = 
-20.63) and is indicated by an arrow in Fig. 8 and the difference between the distributions is 
not at all significant, as confirmed by the KS test which gives a probability of p=0.77 that 
the two distributions are similar. However, if we eliminate 5 ellipticals (the 4 brightest and 
the faintest one) in the high density sample and the 2 faintest ellipticals in the low density 
one, so as to make both samples to cover exactly the same range in luminosity, we find a 
= 3.7, cr(a Q )= 32.1 km s" 1 for the unweighted fit and a = 3.9, cr((To) = 32.2 km s _1 , for the 
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weighted fit in the high density sample, a = 2.7, er(<xo)= 69.9 km s _1 for the unweighted fit 
and a = 3.6, a(ao)= 68.6 km s _1 for the weighted fit in the low density sample, confirming 
what we have obtained on the whole luminosity range. 



Finally to ch eck how so 



(DR 8) database 



York et al. 



i d can b e considered our resu lt we have inspected the SDSS-III 



(120001 ); lAihara et al.l (120111 ). looking for ctq measures for 



ellipticals in the low and high density environments. Unfortunately those data are available 
only for 14 ellipticals in the high density environment and for 12 ellipticals in the low 
density one and are reported in Tables 5 and 6, in which we list for each elliptical in each 
sample, identificator (column 1), o"o value from the SDDS (when available) with related 
uncertainty (column 2) and the difference (Ado) between SDSS and Hyperleda value for cr 
(column 3). Inspection of data in column 3 reveals a Ado which is in general small, but 
almost always negative in the high density environment (<o"o> = -13.2 km s _1 , <o"o>rms = 
9.8 km s _1 ) and that is larger and more spread around the zero in low density environment 
(<a > = 5.9 km s _1 , <a >R M s = 24.2 km s _1 ). 

Replacing SDSS a (when available) to Hyperleda ones gives L B oc cr^ 8 for the high 
density environment (in both weighted and unweighted fit) and L B oc ctq' 3 , L b oc <7q 4 for the 
low density environment (respectively for the weighted and and unweighted fit). Average 
scatter of the data gets somewhat worst, especially in high density environments. We find 
cr(<7o) = 36.5 km s _1 for the high density enironments and cr(o"o) = 68.1 km s _1 for the low 
density ones (to be compared with 33.7 km s" 1 and 67.0 km s _1 which are the corresponding 
values obtained for the weighted fits when using only data from Hyperleda). Despite of the 
slight increase of cr(o"o) which is likely to be due to the non homogeneity between the two 
distinct sets of data (as proved by the shift in o"o reported in Tables 5 and 6, column 3), the 
environment effect on o~(o~o) is confirmed. 

Table 7 summarizes all the results described above, in column 1 we indicate the 
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sample kind (overdense and underdense stand for high density and low density) and within 
parenthesis its caractheristics related to a (either small error or data from the SDSS), in 
column 2 the sample size, in column 3 the kind of fit (either unweighted or weighted), in 
column 4 the FJ relation slope and related uncertainty, in column 5 the average dispersion 
of the data around the fit. 

Finally, in analogy with what we have done for the whole sample (see sect. 3, Fig. 5 
and Table 4) we have checked for possible maximum effects due to me errors on the FJ 
relation in the high and low density environment. The procedure is exactly the same but in 
this case we have generated only 30 simulated samples by random addiction or subtraction 
of the real error on m-Q. The results are shown in Fig. 9 showing the normalized distribution 
of (t(o" ), a and Act for the high density (upper panels) and low density (lower panels) 
simulated samples. The arrow on each plot indicates the value obtained in the real case. 
From Fig 9 we see that errors on luminosity would produce a general degradation of the fit 
quality (particularly evident in the possible large increase of a and Aa for the low density 
sample), but that however the difference in a(o"o) would be mantained. 

6. How much overdense and underdense have to be the environments? 

In the previous section we have shown that ellipticals in high density environments 
display a significant reduction of the FJ scatter, when compared to ellipticals in low density 
ones. Both kind of environments have been selected objectively requiring a number of 
neighbours (N ne i g h) equal or larger or equal or smaller than 5 times or 0.2 times the median 
value of the N neig h distribution (computed separately for each distance bin). 

We now reduce and increase progressively the overdensity and underdensity value (i.e. 
the multiplicative factor that we have applied to the median value of N neigh ) to check the 
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level of densities at which the difference in the FJ scatter holds. 

The results of this test are shown in Tables 8 and 9 where for each value of the 
overdensity or underdensity (column 1) we indicate the number of ellipticals in each sample 
(column 2), the kind of fit (column 3), the slope with its uncertainty (column 4) and the 
average scatter a(a ) of the FJ relation with related uncertainty(column 5). Tables 8 and 9 
allow one to follow the increase of cr(<Jo), as overdense environment becomes less and less 
dense and, complementary, the decrease of cr(<7o), as underdense environment gets more 
and more dense. From Table 8 we see that a (cro) mantains its small value for overdensities 
down to a value of 3.5, that it is still small, even if somewhat increased, when overdensity 
is equal to 3 and that then it starts to increase more rapidly to reach the characteristic 
value displayed by the whole sample at an overdensity of 1.5. From Table 9 instead we see 
that the dispersion is already below 60 km s" 1 at an underdensity factor of 0.25 and that 
it decreases progressively remaining just above the a(a ) of the whole sample when the 
underdensity factor is equal to 0.75. 

This progressive increase/decrease of a (er ) with decreasing/increasing density in 
overdense/underdense environments gives more strength to our result confirming an effect 
relating environment to the FJ relation scatter. 

7. Conclusions 

Using a sample of 384 nearby elliptical galaxies and objectively estimating their 
environment on the basis of the number of neighbours within the typical galaxy cluster and 
group scale we have provided evidence for an effect relating the intrinsic scatter of the FJ 
relation to the environment. We have shown that the scatter of the FJ is reduced to almost 
half of its value when ellipticals in highest overdensities are compared to ellipticals in 



-18 - 



less-density environments, that the effect is not induced by luminosity differences between 
the samples and that it holds for overdensities ranging between 3.5 and 5 the median value 
of the number of neighbours distribution. Besides indicating a rather simple and quite 
natural way to reduce the large scatter affecting the FJ relation, our result, if confirmed 
on larger samples, is very likely to open an interesting perspective for models of galaxy 
formation. 
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Bin 


A v v 
km s -1 


Nell. 


<M B > 


med(M B ) 


I 


2000 <v r < 4000 


61 


-20.25 


-20.21 


II 


4000 <v r < 6000 


136 


-20.59 


-20.67 


III 


6000 <v r < 8000 


120 


-20.98 


-21.02 


IV 


8000 <v T < 10000 


67 


-21.39 


-21.41 



Table 2: The effect of distance on luminosity on our sample. 
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Sample 


* en. 


fit 


a 


(j((Jn) 

km s 1 


whole 


384 


unweighted 


o.l 


51.9 ±1.2 


bright 


206 


unweighted 


4.7+2-! 

1 —0.6 


49.5 ± 1.8 


faint 


178 


unweighted 


2 7 +0 n i 

1 -0.3 


52.9 ± 1.7 


whole 


384 


weighted 


4.018:1 


51.1 ± 1.2 


bright 


206 


weighted 


5.6+°;? 


49.3 ± 1.8 


faint 


178 


weighted 


q n+0.5 
°- z -0.4 


52.3 ± 1.7 


whole (Ao-q/o-q < 0.1) 


314 


unweighted 


q O+0.2 


48.4 ± 1.0 


bright (Ao-q/cto < 0.1) 


174 


unweighted 


5 4+i-o 

°-^-0.8 


47.9 ± 1.5 


faint (Aa /a < 0.1) 


140 


unweighted 


3 +0 - 5 


47.2 ± 1.1 


whole (Ao-q/o-q < 0.1) 


314 


weighted 


4 1+0-3 
^•-■--0.2 


48.1 ± 1.0 


bright (Ao-q/o-q < 0.1) 


174 


weighted 


5-8lag 


48.0 ± 1.5 


faint (Ao-q/o-q < 0.1) 


140 


weighted 


o q+0.6 
O-O-0.4 


47.0 ± 1.1 



Table 3: FJ relation parameters for ellipticals in the whole sample and in the bright 
(Log(L B /L B0 ) > 10.5) and faint ( Log(L B /L BQ ) < 10.5 ) subsamples. 
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Sample 


ArriB 


<a{a )> 
km s -1 


(J ( 'o)rms 
km s -1 


< a > 


OiRMS 


< Aa> 




whole 


0.25 


52.6 


0.6 


4.3 


0.1 


0.52 


0.03 


whole 


[0.03 - 0.72] 


52.8 


0.7 


4.3 


0.1 


0.52 


0.03 


bright 


0.25 


49.8 


1.2 


5.9 


0.5 


1.7 


0.3 


bright 


[0.03 - 0.72] 


50.1 


1.3 


5.7 


0.5 


1.6 


0.3 


faint 


0.25 


54.5 


1.5 


3.5 


0.3 


1.1 


0.2 


faint 


[0.03 - 0.72] 


54.2 


1.8 


3.7 


0.3 


1.2 


0.3 



Table 4: Average values for the FJ relation parameters derived accounting for possible errors 
on tob- 
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Identificator 


<r (SDSS) 


Aa 




km s _1 


km s _1 


PGC 6945 


- 


- 


NGC 704 


- 


- 


NGC 3837 


- 


- 


NGC 3842 


291 ±5 


-24 


NGC 3862 


260 ±5 


-11 


NGC 4261 


- 


- 


NGC 4473 


- 


- 


NGC 4816 


- 


- 


PGC 44137 


- 


- 


NGC 4839 


269 ±5 


-16 


NGC 4842A 


208 ±4 


-8 


PGC 44367 


144 ±3 


-12 


PGC 44467 


- 


- 


NGC 4860 


263 ±4 


-13 


IC 3959 


205 ±4 


-3 


NGC 4864 


- 


- 


NGC 4869 


_ 


_ 


NGC 4874 






NGC 4881 


193 ±3 


-7 


NGC 4882 


149 ± 3 


-14 


NGC 4884 






NGC 4906 


173 ±3 


1 



Table 5: Ellipticals in the overdense environment 
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Identificator 


<r (SDSS) 


Aa 




km s _1 


km s _1 


IC 4041 


119 ±3 


-17 


IC 4045 


208 ±3 


-9 


PGC 44848 


173 ±3 


-40 


NGC 4926 


264 ±4 


-12 



Table 5: continued. 
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Identificator 


<r (SDSS) 


Aa 




km s _1 


km s _1 


NGC 631 


- 


- 


NGC 810 


292 ±5 


34 


NGC 1226 


- 


- 


UGC 3549 


- 


- 


UGC 3844 


189 ±3 


39 


NGC 2474 


- 


- 


NGC 2800 


- 


- 


NGC 2954 


188 ±3 


-28 


UGC 5313 


63 ±4 


-11 


IC 590 


296 ±5 


23 


NGC 3392 


165 ±2 


6 


NGC 3731 


154 ±3 


-19 


NGC 4187 


277 ±5 


-23 


NGC 4272 


- 


- 


UGC 7813 


256 ±4 


-16 


NGC 5583 


- 


- 


NGC 5628 


254 ±4 


29 


NGC 5771 


122 ±2 


22 


IC 1101 






NGC 6020 


205 ±3 


15 


NGC 6051 






IC 1211 







Table 6: Ellipticals in the underdense environment 
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Table 6: continued 



Identificator cr (SDSS) Aa 
km s _1 km s _1 

NGC 6442 
NGC 6515 
NGC 6575 
NGC 6697 
NGC 6702 
IC 1317 
NGC 7052 
NGC 7360 
NGC 7512 
PGC 71599 
NGC 7735 
NGC 7751 
NGC 7785 
NGC 7786 
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Sample 


N ell . 


fit 


a 


o"(o"o) 
km s _1 


overdense 


26 


unweighted 


O O+0.5 

O.8_o.4 


33.7 ± 1.4 


overdense 


26 


weighted 


q q+0.5 


33.7 ± 1.4 


overdense (with 14 SDSS cr ) 


26 


unweighted 




36.6 ± 0.9 


overdense (with 14 SDSS a ) 


26 


weigthted 




36.5 ± 0.9 


underdense 


36 


unweighted 


q q+1.4 
°-°-0.7 


67.7 ± 3.7 


underdense 


36 


weighted 


3-811:? 


67.0 ± 3.7 


underdense (Aa /a < 0.15) 


30 


unweighted 




68.5 ± 2.9 


underdense (A<7 /<7 < 0.15) 


30 


weighted 


3.4115 


67.9 ± 2.9 


underdense (with 12 SDSS a ) 


36 


unweighted 


o 4 +i.5 


69.0 ± 2.7 


underdense (with 12 SDSS cr ) 


36 


weighted 


4 0+2.1 


68.1 ± 2.6 



Table 7: The FJ relation in overdense and underdense environments. 
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overdensity 


Well 


fit 


a 


o"(o"o) 
km s -1 


4.5 


37 


unweighted 


3-6^ 


32.0 ± 1.4 


4.5 


37 


weighted 


q 7+0.4 
°- ' -0.3 


32.0 ± 1.3 


4.0 


45 


unweighted 


q 7+0.4 
°- ' -0.3 


34.1 ± 1.2 


4.0 


45 


weighted 


q O+0.5 

<J-°-0.3 


34.0 ± 1.2 


3.5 


56 


unweighted 


3 6 +0 - 3 


33.2 ± 1.1 


3.5 


56 


weighted 


q O+0.4 


33.0 ± 1.1 


3.0 


73 


unweighted 


3 6 +0 - 3 

°-°-0.3 


35.6 ± 1.2 


3.0 


73 


weighted 


q O+0.4 


35.0 ± 1.2 


2.5 


88 


unweighted 


q O+0.3 


41.1 ± 1.1 


2.5 


88 


weighted 


4 +0 - 4 

^• u -0.3 


40.7 ± 1.1 


2.0 


108 


unweighted 


o q+o.4 

°- y -0.3 


45.6 ± 1.9 


2.0 


108 


weighted 


4 1+0-3 


45.2 ± 1.9 


1.5 


138 


unweighted 


q O+0.4 


51.4 ± 1.8 


1.5 


138 


weighted 




50.9 ± 1.8 


1.0 


198 


unweighted 


3.61°;! 


51.1 ± 1.5 


1.0 


198 


weighted 


4-11^ 


50.2 ± 1.5 



Table 8: Relaxing the overdensity. 
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under density N c u fit a cr(co) 

km s _1 

0.25 78 unweighted 3.4±°;£ 59.3 ± 2.6 

0.25 78 weighted 3.9±g;? 58.8 ± 2.6 

0.40 99 unweighted 3.5+^ 58.1 ± 2.4 

0.40 99 weighted 3.9ljjj 57.6 ± 2.4 

0.50 127 unweighted 3.ll^ 55.9 ± 2.6 

0.50 127 weighted 3.6^ 55.2 ± 2.5 

0.75 158 unweighted 2.9±g;l 53.8± 2.2 

0.75 158 weighted 3.4+°;^ 53.2 ± 2.2 



Table 9: Relaxing the underdensity. 
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Fig. I. — The distribution of the whole sample ellipticals (filled circles) in the v r , M B plane. 
All data lay above the curve representing the increase of the minimum observable luminosity 
(Mb) with distance for galaxies having m B < 15.5. (This latter is the limit that we have 
imposed to the sample). 

Fig. 2. — Relative contribution to the absolute blue magnitude (Mb) distribution of the 
whole sample (dotted) from ellipticals belonging to 4 bins of increasing distance (see Table 
2). 

Fig. 3. — The FJ relation for the whole sample. The green continuos line is the best fit for 
the whole sample (Lb oc (Tq°). The red dashed ( Lb oc <Jq 6 ) and the blue dotted (Lb oc o"q' 2 ) 
lines are the best fits for the low ( Log(L B /L B 0) < 10.5 )and high (Log(L B /L B 0) > 10.5) 
luminosity ellipticals. All fits have been derived accounting for errors on er . 

Fig. 4. — The distribution of relative errors on er (Aa /a ) for ellipticals laying above 
(dotted) and below (continous) the best fit line in the whole sample (upper panels) and 
in the high and low luminosity subsamples (middle and lower panels). Left panels refer to 
unweighted fits, right panels to fits weighted for errors on er . 

Fig. 5. — The effect of the error on itlb on the FJ relation. Upper panels refer to the whole 
sample, middle and lower panels to the high and low luminosity subsamples. Each distri- 
bution (dotted or continuous histograms) refers and is normalized to a set of 300 simulated 
samples that we have obtained randomly adding or subtracting to each m-Q its true (dot- 
ted histogram) or average (continuos histogram) error. The arrow on each panel indicates 
corresponding values obtained for the real sample and subsamples. 

Fig. 6. — The number of neighbours (iV nei g h ) distribution for ellipticals in each of the 4 
distance bins in which we have divided the whole sample (cfr. Table 2). 
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Fig. 7. — The FJ relation for ellipticals in high density (red filled circles) and low density 
(blue open squares) environments. Overimposed are the corresponding best fit lines (red 
continous L B oc o"q' 79 and blue dotted L B oc ctq' 75 ). 

Fig. 8. — Mb distribution of ellipticals in low and high (shaded) density environments. Both 
distributions are normalized to the total number of ellipticals in each sample. The arrow 
indicates the median value of M B which is the same for both samples. 

Fig. 9. — The effect of the error on m B on the FJ relation for the high density (upper panels) 
and low density (lower panels) environment. Each distribution refers and is normalized to a 
set of 30 simulated samples that we have obtained randomly adding or subtracting to each 
tub its true error. The arrow on each panel indicates corresponding values obtained for the 
real samples. 
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